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Abstract 

Entangled continuous variable (CV) Gaussian states with different wavelengths plays a central role in 
recent CV-based approaches to quantum network, quantum information processing and quantum metrol¬ 
ogy. Typically, experiments demonstrating CV entanglement exploit the optical parametric frequency down 
conversion. Due to the probabilistic nature of photon pair generation, the entanglement involving the post¬ 
selection of photonic qubits is limited to at most three colors. Here We theoretically present a scheme for 
the deterministic generation of entanglement among bright multicolor CV Gaussian states from an optome¬ 
chanical system using existing experimental technologies. In our scheme an optical frequency comb is input 
into an optomechanical resonator and then the amplified optomechanical coupling makes multipartite en¬ 
tanglement among them. Our scheme overcomes the limitation of usable frequency of entangled CVs in 
frequency conversion process. It can be extended to generate multipartite entanglement between orthog¬ 
onal modes with a single frequency or between microwave and optical CVs, or even among a microwave 
frequency comb. This optomechanical device can be integrated on a chip. 
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Multipartite entanglement of continuous variable (CV) cluster states is not only of fundamental scientific 
interest [|TI|, but also the key ingredient for quantum information technologies such as universal quantum com¬ 
putation [l2l|3l|4l|, quantum metrology [|5l[^, gravitational wave detection Q, and even quantum network of 
clocks [[51. 

A number of different techniques for the generation of entanglement in the CV regime have been proposed 
and experimentally realized. Entangled CV cluster states are created by combining squeezed Gaussian states 
generated from optical parametric oscillator (OPO) on beamsplitter but is limited to one frequency HOl HH 
dll . Signal and idler photons generated from OPOs during the frequency down conversion (FDC) are naturally 
entangled but limited to two color [[l3l[l4l|T5l[I^. This FDC-based technique is only extended up to entangle 
three-color CVs [[Il[l2l due to two limitations: the probabilistic nature of FDC and the available frequency 
reduced by half in each FDC. We note that Four-wave mixing technique has been demonstrated to be able to 
create quantum correlation among three-color CVs but entanglement is not clear liTSl . 

To date, entangled CV cluster states has only been generated up to two colors in frequency domain via the 
probabilistic FDC recently IfTSlfTSlfT^ . or in spatial modes IfT^ . or in time domain via optical group delay [|9l. 
However, the wavelength or color of entangled CV states is limited by available optical nonlinear crystals or 
materials. In this work, we propose a method using a multimode optomechanical system to create entangle¬ 
ment among multicolor CV cluster state. We apply an optical frequency comb (OFC) of coherent laser field 
to drive the cavity modes on-resonance. In our configuration, the reflected CV Gaussian states off cavity with 
at least ten colors around each cavity mode frequency are entangled. Our scheme overcomes the fundamental 
limitation of mode number of entangled CVs in frequency domain. 


Results 

Theoretical description. The setup for generation of multicolor CV entanglement is illustrated schematically 
in Fig. [Ja). The optical frequency comb is filtered and then incidents to the first optical grating. The spectra 
decomposited by this grating is modulated by a liquid crystal light modulator (FCFM) with fast response 
to the amplitude and phase of each spectral line IIT3l . and then is combined by the second optical grating. 
The modulated optical frequency comb incidents into the optomechanical resonator supporting optical cavity 
multimode, which has resonance frequency coj and intrinsic loss rate of k^j for the jth mode. The jth CV mode 
of optical frequency comb at frequency u>l,; drives the corresponding jth cavity mode with an amplitude of ej, 
and an external coupling rate of /Ce,;- One of two mirrors of optomechanical resonator oscillates with frequency 
Dmo nnd its motion decays with a rate of The mechanical motion of this movable mirror couples to the jth 
cavity mode with a rate of gom,i- The optical fields reflected off the OMR is isolated from the input fields by a 
highly reflective beam splitter. Using this setup, we can create entanglement among the multicolor CV modes 
of OFC. 

Under the rotating-wave approximation of the driving, the Hamiltonian describing the interaction between 
the optical modes and the mechanical motion is liT9l : 

H = Qmob'^b -I- ^ A'.a]a/ -i- ^ gom,ja\a0 + b) 

_^ ■’ ( 1 ) 

-H i ^ - e)aj) , 


with the detuning A'. = ojj - a>i^j. Here aj and b are the annihilation operators of the jth cavity mode and 
the mechanical motion, respectively, gomj is the zero-point optomechanical coupling rate of the jth cavity 
mode. Driven by multi CV modes in OFC, the cavity modes and the mechanical motion can reach the steady 
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Figure 1: (a) Schematic setup for generation of multicolor entanglement in optomechanical systems. The 
spectrum decomposited by the fist optical grating from an optical frequency comb is modulated by a liquid 
crystal light modulator (LCLM) and then is combined by the second optical grating. This modulated frequency 
comb is input into a multimode optomechanical resonator and subsequently is reflected off the optomechanical 
system. The optomechanical resonator creates entanglement between the spectral lines (with frequency ,) 
of the input optical frequency comb, (b) Entanglement of two modes when G = O.SQmo- The plot shows FJ 2 / 2 . 
Other parameters are A = 0, Omo = O-l/c, Qm = 10®, fith = 10^, Ki = 0. 


-——^———. By substituting aj aj + aj, b ^ /3 + b and the output aout,; ^ cioutj + cioutj, we linearize 

the operators to a weak fluctuation displaced by a strong coherent shift aj, (5 and dout,; ifTOllSOll . respectively. 
Normally, IdoutjP, he- the photon number of output field, is much larger than unity, implying a bright output 


beam. Correspondingly, the fluctuation in the X and Y quadratures are AXout,Xm) 
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from environment becomes 


m- After linearization, the Hamiltonian describing the dissipation and noise 


H' = Q^ob^b + Aja]aj + + G*aj)(b^ + b) , 


( 2 ) 


with Gj = gomjbj. Obviously, one can tune Gj to be real by controlling the phase of driving, sj. Throughout 
the investigation below, we assume all the amplified coupling rates Gj real for simplicity. We now study the 
quantum fluctuation in the output using the quantum Langevin equations (QLEs) [I22ll23ll : 

Q-j — {iAj Kj^Qj iGj(b Z?) “1“ 'sj 2/c”e,y^iin,e,y(t) “t" '^2/(’i yHjn.i.yCO , (2a) 

b = -{iO-niQ + ymo)b - i 'Yj + <^ 7 ) + V27m0^in(0 , (3b) 

j 


where the optical vacuum noise entering the cavity through the external optical coupling and the intrinsic optical 
loss channels are included by the random fluctuating inputs ain,e,y(0 and respectively. bin(t) describes 

the mechanical noise applying to the mechanical resonator, y^o is the decay rate of mechanical motion and 
therefore gives mechanical quality factor 2„, = Clmo/y,„o. To solve the QLEs in frequency domain, we define the 
Fourier and inverse-Fourier transformations A(m) =dcoe A{t) and A{t) = dcoe A(co) [| 21 [| 24 [ 

I 25 II . Thus we have the correlation for noises {ain,e,j(co)alne + m'), = 6ji6(a) + 

to'), where 6ji {6(a) + to')) is the discrete Kronecker (Dirac) delta function. Here we assume that the occupancy 
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of cavity modes due to the thermal environment is negligible, = 0. But the optieal 

vaeuum noises from external eoupling is uneorrelated to that from the intrinsic loss, i.e (ain,ej(ru)a[j^ ■/(cu')) = 0 

for any j and 1. We also have {bin(a))bl^(a)') = (ilth + l)d(m + to') and {bl^(a>)bin(co') = nthd(m + to'), where iith = 
(^^nn„,olKBT _ thermal occupaney of mechanical resonator at temperature T. But {binito)bin{to') = 0 

and {bl^(to)bl^{to') = 0. The solution to the QLEs, Eq. in the Eourier domain is given by [I2n[i4ll25]l 


aj(to) = 
b(aj) = 


^j2K^jXR,j{to)am,^j{to) + ^j2ixjXR,i{.oo)am,tj{to) - iGjXRjioj) [b\to) + b(a>)^ , 
^f2^bin(oJ) - i Gj [a^jito) + a/m)] 
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where the optieal eavity suseeptibility for the jth mode is Xrj((^) = 1/[*(A/ - to) + Kj]. The mechanical 
susceptibility Xm((^) conneeting the meehanieal response to the environmental noise and the optical vacuum 
fluctuations incident on the optieal eavity takes the form;^f“/m) = f(f^mo-<^)+7mo + Zj Gj XR,j(<^) “T«y(“^)]- 
The fluetuations in the output fields can be ealeulated from the input-output relation [|26l 123 aout./oj) = 
-a-mfijioo) + yl^eQjioo). Then we have the output noise eorrelations [fTl l2TH28]l 
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“1“ 2 'y„j 0 ^KgjKglGjGl^j{oo)^i (to) j^l 7 m(aj)| (f^m “1“ 1) “1“ \Xin( 0^)1 
+ 2 yfl^,GiGj% [(2kjXrj(oj) - l) F:xoo)s;(aj)xlj(oj)e-'^] , 

where = Xm(-^) -Xm(oj), ^j(oo) = e'V*/-m) - e-'^XRj(^) and Ej(to) = e'\\.{-oo) + e-'^XRj(<^)- 
means the real part of number x. 

One of the eriteria for analyzing multipartite entanglement of CV modes is Duan eriteria Il29ll and its 
extension [1^ written direetly in terms of these eorrelations, as sums of varianees Vji between modes j and /: 

Vf\to) =<(AXy(m) - AAKm)) (AXy(-m) - AX^-m))) 

+ <(Ayy(m) + Ay;(m)) (Ay/-m) + AyK-m))) > 2 , 


for j 1. It suffiees to demonstrate M-partite entanglement as long as the above inequality is violated for j < I 
and j,l 6 {1,2, • ■ • ,M] (U |3T1. The degree of violation indicates the degree of entanglement. Note that the 
inequality Eq. I^is not necessary optimal for searehing the largest entanglement but its violation is a sufficient 
criterion for inseparability. 

Next we present a simply closed formula for estimation of the eorrelation and variance for the ease all 
eavity modes are identical, i.e. Ay = A, Gy = G, /Ce.y = Ks, and = /Ci. Although bipartite squeezed states have 
been proposed in optomechanical systems by applying the red-detuned and blue-detuned multitone driving 
simultaneously Il3^ . it is hard to be scaled for multipartition entanglement. Here we take the scheme using a 
single zero or small detuned driving recently demonstrated in experiments for squeezing a single optieal mode 
in optomeehanieal systems [121112^[^ . Therefore, we take A « 0. We define 6+ = Qmo ±aj and 6^^ = 6Z^ +d/, 
and assume that |d| » ymo- Under this eondition, for eavity modes with identical decay and optomechanical 
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coupling rates, we have 
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with 77 = Ke/K and fn 


4G^ 


According to Eq. 


if one quadrature, e.g. Xout, is correlated, another 
ig. [^b), the minimal variance is obtained around 


quadrature, Tout, will be anticorrelated. As shown in 
\6\ = nl2, 3;r/2. However the frequency width of entanglement is small (also see supplementary material). In 
contrast, the optomechanical system has a broad width where entanglement is realized when \6 + ;r/ 2 | = pi 14. 
Next we focus on \6\ = n/A. In the region of |m| ~ such that 6^{52^ + 6~^) « I, the variance as sums of 
correlations is 
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The first term is the Duan’s bound. The entanglement between the jth and Ith modes requires |d| > Mfnieas/2 + 


Q,n 

MY 


, but the minimal variance, 
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meas + AYI^q^. Interestingly, the largest achievable entanglement is independent of the optomechanical 


4rr meas +400,0 ^ 


> 2 - 77 /M, is available at the optimal frequency |d, 


optl 


coupling strength, but is limited by the number of involved cavity modes. 

According to our analysis above, the perfect coupling regime [l2n . i.e. 77 = 1, is preferable for our aim to 
create entanglement. In this regime, the reflectivity of the coherent driving fields are almost unitary. There¬ 
fore, the reflected beams include many photons implying bright output beams. The violation of inequality 
Eq. (30) implies the entanglement of strong CV Gaussian states. Although one-mode squeezing of light has 
been demonstrated in optomechanical systems in both the unresolved-sideband regime lilTl and the resolved- 
sideband regime Il28l . we will focus on the unresolved-sideband regime and the perfect coupling regime in our 
discussion below. We also take values for parameters, A = 0, = 0.1/c, Qm = 10^, fith = 10^, Ki = 0. 


Entanglement of two CV modes. When two cavity modes are identical in detuning. A, = A = 0, optomechan¬ 
ical coupling strength, G, = G, and decay rates, /Cej = /Ce and = Ki, the variances Vji are equal for any j and 
1. If we drive the cavity modes properly, we are able to create entanglement of multicolor CVs. Eor example, 
the variance can be below the Duan’s bound (blue region in Eig. [TJb)) in the centre of the reflected CV modes 
corresponding to the driving modes in the output, leu - < 0.95^2^ for 6 = -n/4, -3nlA. Entanglement also 

can be obtained over the frequency region of |eu - > 1.03Gm for 6 = n/4,3n/4. The minimal variance can 

be down to almost 1 at 6 = nl2, implying a violation of Duan’s bound by 3 dB. 


Four-color entanglement. We now examine how the random varying optomechanical coupling and decay 
rates, which are hard to control precisely in experiments, affect the degree of entanglement. In Eig. we 
calculated all variance Vji for j < I for four output CV modes. Pink lines show the variance for identical pa¬ 
rameters and samples using Eq. (|^ in Eq. ( [SO] ). To analyzing the influence of random variation of parameters 
in various samples, we present the variances of ten samples with independent normally distributed random 
optomechanical coupling rate, G, and decay rate, k^j with variance of 10% of their mean value. It is found that 
all six variances (7 < 1) can be as low as 1.75, yielding a considerable four-color entanglement. 
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Figure 2: Four-mode entanglement evaluated over ten random samples. Parameters {Gj and Kej) of samples 
vary within 10% of mean value with a normal distribution. Pink line indieates the varianee for samples with 
identieal eavity modes. Dashed red (green) line shows the Duan bound [|29l [301 (the limit given by Vniin)- 
G = O.SDmo^^ = ^/4. 


Ten-color entanglement. Ten-eolor entangled CV eluster state is illustrated in Fig. Figure [^a) shows the 
variance evaluated by Eq. (30) using Eq. ([^ between modes 1 and 2 among ten CV modes. The pink line 


shows the variance for identical cavity modes. In the case of identical cavity modes all variance VY 2 with 


(iO 


j I are equal. The blue spots show the variance of ten random samples with random parameters Gy 
and /Ce,y taking a normal distribution with variance of 5% of mean value. Obviously, this amount of variance in 
parameters only changes the entanglement nature slightly. The variance between any two of ten cavity modes 
in an optomechanical resonator with identical Gy and /Ce,y for one sample is shown in Eig. [^b). The minimal 
variance Vmm is dependent on the optomechanical coupling, G. It decreases from 2 to a limit of 2 - 77 /M as 
G increases, see blue line. Considering that the characteristic of solid-state devices is diverse from sample 
to sample, we examine the variance of each CV mode pairs over one hundred random samples, of which the 
parameters have the same mean value, but is distributed normally with variances of 5% (see Eig. |^c)). Clearly, 
all variances, are below the Duan bound, although they may vary over a small range. The average is about 
1.9. It implies a ten-color entangled CV cluster state. 


Experimental implementation. We take a Eabry-Perot cavity containing a Si 3 N 4 membrane in the middle 
as our optomechanical implementation. The whole system is at low temperature, T ~ 500mK [|34l|351. Eor 
simplicity we assume that this EP cavity supports 10 modes with different resonant frequencies but identical 
decay rates Kjln = IMHz [[34l and Kijln negligible [l34l . The membrane we use has 2m = 4 x 10® and 
Dmo/27r = 134kHz Il35l [3^ corresponding to the thermal phonon number of iith = 8 x 10^. We drive the 
cavity modes with a frequency comb from on-chip optomechanical resonators OTll^ . Under on-resonance 
driving at each mode with input powers of Pin = 80 qW exciting the photon number of ~ 1.4 x 10* in 
each cavity mode around A ~ 1064nm. Such driving is strong enough to amplify optomechanical cou¬ 
plings to O.lDmO- In this optomechanical system the minimal variances of ~ 1.89 < 2 is obtained 
at to - tOLj ~ 1.052Dnio U 6 {1,2, • • • , 10}). Therefore, we create an entangled ten-color CV cluster state. 
If we apply our scheme to the microwave-optomechanical system developed by Andrews et al. [|39l . we 
can create entanglement between microwave and optical signals by reducing the variance to ~ 1.3 
(Qmo/27r = 380kHz, 2^ > 10^,Dnio//<^ ~ 0.25, T = 40mK, ilth = 2.2 x 10* and G = O.lDmo)- 


Discussion 


6 



























Figure 3: Entangled ten-mode CV eluster states, (a) Varianee of Parameters of ten samples, Gj and 
Ksj, take a normal distribution with varianee of 5% of mean value, (b) Varianee of (j + 1) (pink 
lines) of one sample with identieal parameters, Gy and k^j, as a funetion of to for various eoupling rates 
G = {0.02,0.05,0.08,0.1,0.2,0.3,0.4} x Blue line shows the positions of smallest varianees. (e) Varianee 
matrix, for 7 , Z 6 {1,2, • • • , 10} at tu = 1.04QniO- Bars show the minimal and maximal varianees evaluated 
over one hundred (5 = 100) samples with random coupling rate, G^j, and external coupling rate, varying 
within a normal distribution with variance of cr = 5%. G = O.lflmo, ^ = -Tt/A. 

In summary, we studied the multicolor entanglement of bright CV states from optomechanical systems. En¬ 
tanglement up to ten CV modes has been demonstrated, while the output fields have many photons. It is found 
that the largest violation of Duan’s bound among two-color CV modes is bounded by a limit of 3dB (|0| ~ 7 r/ 2 ). 
To create multicolor entanglement our on-chip optomechanical setup only requires on-resonance driving of the 
corresponding multimode of cavity. Our proposal does not rely on the probabilistic wave mixing process and 
therefore can also be extended to generate entanglement between microwave and optical photons [l39l l40ll or 
even in microwave frequency comb [I4ni4^ . Our work provides routes towards microwave or optical quantum 
frequency comb. It may enable quantum frequency comb-based applications in quantum information process¬ 
ing and quantum network. 
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Supplementary material: Deterministic generation of bright multicolor entanglement from 

optomechanical systems 


In this supplementary material we present the detailed derivation of the output noise power spectral density 
of continuous variables (CVs) and the formula for their variance between two different cavity modes. 


1 Model 

1.1 Hamiltonian and quantum Langevin equations 

Our setup for multicolor entanglement of CVs is schematically illustrated in Figure 1(a). An optical frequency 
comb (OFC) laser incidents on the first optical grating and then is decomposited into a set of spectral lines 
with carrier frequency ujlj. Each spectral line is modulated in amplitude and phase by the liquid crystal light 
modulator for proper driving of cavity modes. These spectrally modulated OFC is combined by the second 
optical grating and then is applied to drive the cavity modes with rates /Ce^ . The OFC field is reflected off the 
cavity to the output. The quantum noise in the output fields is squeezed due to the optomechanical interaction. 
The Hamiltonian governing the evolution of multimode optomechanical system in Fig. 1 takes 

H = Q.^Qb^b + ^ A'.aja, + ^ gom,ja]aj{b^ + (?) + z ^ - s*aj) , (9) 

,/■ i j 

where A' = coj - coi^j is the detuning between the jth cavity mode and its driving, gomj is the optomechanical 
coupling between the cavity mode a, and the mechanical motion b. 

The quantum Langevin equation reads |23l 

hj — ^gom,jib + fi) + 'sJ'l.Kg jSj + 

b — {iClffiQ + yinO^b i ^ g^mjCljO-j + mobm (0 5 ^ ^ 

j 

Where and y^o are the intrinsic loss rate of the j’th cavity mode and the mechanical decay rate of mechanical 
resonator. All loss rates are necessarily accompanied by random fluctuating inputs ain,ej, tZin.i,;, and bin{t), for 
optical quantum noise from the external coupling channel, intrinsic loss channel and mechanical noise. We 
have <ain,ej) = = 0 and <(?!„) = 0. 

In the steady state, we have 


fi = {b)ss = 


i gomjl^jl 

idfnO TmO 


-y 

XnO ; 


/ J gom,j\^j\ 5 


iAj + Kj 


(lla) 

(llb) 


with Aj = A'j + 2gom,jj3. 

Now we linearize the operators as aj aj + aj and b ^ /3 + b lfT9ll20ll . and then have the Hamiltonian after 
the linearization reads as 

H' = Q^ob^b + + G*aj)(b^ + b) , (12) 

,/ 
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with Gj = gomjcij- Then the quantum Langevin equations (QLEs) for a and b becomes 

b = -(iQmo + 7mo)b - i Y^iGja'^j + G*aj) + V2rm0^in(0 , 


(13a) 

(13b) 


with the effective detuning becomes Ay = cuy - ojlj + '^gomjP and the total decay rate of cavity mode is 
Kj = Ke j + Kij. Below we assume all Gy are real because we can tune the detuning and the driving sj. 

To solve the Langevin equation in frequency domain, we define the Fourier transformation 


— r 

The optical cavity susceptibility for the yth cavity mode is 


Ait) 

Moj) 


1 r 


d(oe '‘^'A{(jj) , 
dcoe“^*Ait) . 


XRji^o) 


1 


Z(A; - io) + Kj 


The mechanical susceptibility ;ym(m), 

Xm(oj) = KO-mO - m) + y„,o ^'YjG] [xRjico) -XR,ji-<^)\ , 


(14a) 

(14b) 

(15) 

(16) 


connects the mechanical response to the environmental noise and the optical vacuum fluctuations incident on 
the optical cavity. Applying the relation [A(-a»)]^ = A^(m) and [A^(-m)]^ = A(m), we have ll^ 


ajiio) = ^j2K^jXR,jiio)am,^jiio) + ^j2^jXR,iio))am,iji(A) - iGjXRjioj) [b\a)) + h(m)] , 
y/^ymobinioj) - i Z , Gj [aj(m) + ay(m)] 


bio)) = 


(17a) 

(17b) 


z(f^m0 -co) + y^o 

At temperature of T, the correlations for noises are {ain,ej(co)aln e /(^O) = («o,;+l)dyzd(a»+a»'), {am,i,j(oj)(4n i /(^')) 
(hoj + l)dy/d(m + a>'), where 6ji (6(a) + o)')) is the discrete Kronecker (Dirac) delta function and ho,y = 
{e^jlKsT _ Typically, hoj » 0 at light frequency. We also have the correlation for mechanical noise, 
{bin((jo)b\j^(x)') = (hth + l)6(a) + m') and {b\j^oS)b{n(<jo') = nthd(m + oj') with hth = - l)"k no.y and hth 

are the thermal photon and phonon occupancies, respectively. 

By manipulation of Eq. ( [TT] ) the mechanical motion in frequency domain can be expressed as the response 
to mechanical noise and the optical vacuum fluctuations through the optomechanical coupling 


b(co) = yl2y^oXm(co)bin(iO) - iXmiOj) E Gj [xRj((^)amfi,j((^) 

j 

- iXm(oj) E Gj [xRj(a))ain,ij(a)) +x*Rj(-oj)al . j(a))] . 


(18a) 


Substituting Eq. 18 into aj(oS), we have 

ay(m) = ^j2KejXR,j(^)ain,ej((^) + ^j2l^jXRj(<J^)ain,i,j(co) 

- i yj2y^oGjXRj(co) [xm((o)bin(oj) + Xm(-(A)b]^(a))\ 
+ 


GjXRj((^){Xm(-<^) -Xm(C0)] ^ ^'J^n\XRn(0J))ain,e,n(<A) + X*Rn(-^Win,e,n(^^^ 

n 

GjXRj(oA){Xm(-<^) -Xm(co)] ^ V^lTR«(^i>)<3,>,,,>(m) + X*Rn(-^')<^]nJ,n(^^^ • 

n 

Using this equation we are able to calculate the power spectal density and the variance of the output noise. 
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1.2 Output field 

The output field eonneets to the eavity mode field through the input-output relation, aoutjioj) = -«in,e,y(<^) + 
yJlKejajico). Defining F^ico) = - Xm(oj) (Note that = -Fm(<u)), then we have the output field 


^out,y(^) [ 1 gj(cu) -|- '2. "\Jl(^jKijXRj{(F)^in,i,j^^^ 

- i V 2 ^ ^^KsjGjXRjico) [xm(oj)bin(oj) 

Z l sje i* 

Gn J — \XRn{iF)ain,e,n{to) + 


( 20 ) 

-r 2 KejGjXRj{iO)F,n{.CO) ^ Gn J— \XRn{co)ain,i,n{<^) +TL(-^)«L , „(<^)] • 

n V 

To study the entanglement of the output field, we define the fluetuation operators for X and Y quadratures 

asClEIl 


/(^) 


ATout,/(t^) 


e ’^aoutjico) + e'^al^^jico) 


V2 


e ‘^aoutjioj) - .(m) 


i'e^t 




/V2 


(21a) 

(21b) 


1.3 X quadrature 

We define X quadrature of the output as 


A2fout,/(<u) — 




V2 


^ ,, A A , ^ Jd{ 


= -^e ' {2Ke,jXR,jito) - l]ain,e,;(m) -h —e' {2Ke,jX*Rj{-co) - l]a;„^^.(m) 

-r ^2KejKije~'\Rj{co)ainjj{LO) + ^j2KejKijF\lj{-(^))a]^ij{(x)) 

+ i V2/Ce,y V^Gy^y(m) \xm{<^)bin{iO) + Xm(.-<^)b]n{<^)\ 

- yf2KejGjFm{C0Xj{C0) V G„ j—\XRn{<^)ain,e,n{(^) + X*Rn(-^)(^]n e 

- ^KejGjF,„(0jXj(0j)^Gn J^[XRn(Oj)ain,i,n(<^) + TL(-^)«L / „(^)] > 

„ V '^ej 


with ^y(m) = [e‘^x*Rj(-0J) - e and ^y(-m) = -f*(m). 


( 22 ) 
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Thus, the correlation of x quadrature between the jth and /th cavity modes can be evaluated by 


^)) ~ ^YmO 'sj^ej^elGjGl^j(o))^i (oj) (^m 1) Ufm( ^)l 

1 2-1 2- 
+ -12k,j;(:rj(oj) - 1| (n^j + 1)6ji + -\2KejXR,ji-^) - 1| nojdji 

- e~‘^ ^KejKelGiGj \2KejXR,ji<^) “ l] + 1) 

- e'^ ^KejKelGiGj [2K,jXRji-<^) “ l] K(<^)C(<^)XRj(-<^)no,j 

- e‘^ ^jK.jKelGiGj \2KejX*R,l^<^"> ~ l] F,n(oj)^j(co)XR,l(^)(^o,l + 1) 

- -^KejKeiGiGj - 1] F^(a>)^i((i>)xRj(-(i>)ho,i 

+ 2KejKij\xR,j(co)\^inoj + l)d^7 + 2KejKij{)CRj(-aj)fnoj6ji (23) 

- 2K,jy[F;jF-e-'^GiGjWRj(o^^^^ + D 

- 2kij 

- 2Kij^jKejKeje'^GiGj\XR,lico)\^Fm(co)^j(co)ino,l + 1 ) 

- 2a:,-,/ ^Ke,lKeje '^GlGj\XRji-<X))\^Fmi(X))^ji(X))ho,l 

+ 2 ^KejKe^lGiGj\F„{<j))\^^j{(j))Ci{iO) + 1) 

n 

+ 2 ^JKeJKeJGlGj\F„(<X^)\^^j{^0)Cli<^) + ^/,n)G^iTR,n(-^i>)P«o,n • 


Merging some terms, we have 


(AXout,/m)AXout,/(-m)) - -(2noj + l)dy/ 

^YmO "\fFeJFelGjGi^j(oS)^l [^Ufm(^)l (^m 1) \Xmi. ^)l 

- xfJ^iGiGj [2kjXrj(io) - l ] F*Jco)^;(aj)xlj(aj)(n,j + 1 ) 

- ^KejKeiGiGj [2KiXRj(a)) - l] Fm(a))^j(a))xRj(oj)(no,i + 1) 

n 

- e~‘^ ^jKejK.iGiGj \2KjXR j{-ix)) - l] F*Jio)Ci{io)xR,j{-co)no,j 

- e‘^ ^KejKeiGiGj [2 kiXr,i(-co) - 1] Fm(co)^j(co)x*R/-co)noj 

+ 2 ^K,jK,jGiGj\FJa))\^Fj(^XH.^') Xi XnGl\xR,n(-(o)\'^no,n • 
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1.4 Y quadrature 

Similar to X quadrature, we define Y quadrature as 


AYontjico) 

1 




i^^2 


6 l2.Kg j(^CO) 


1 


-e''^[2KejX*R,i(-io) - 1]«L 


L e,/xviv,/v—/ —/ /—'■ 

/V2 

- i ^j2KejKije~'^XRji<^)am,i,jioj) + i ^J2KejKije‘^X*Rj(-^Wm,i/co"^ 

- ^J 2 k^j ^ly^GjSjiio) \xm{(^)bin{(o) + Xm{-<^)b]j.(^)\ 


m,ej' 



{XRn{C0)ain,e,n{^) + g «(^)] 


- i ^/2KejGjF„t(a))Ej(a)) ^ G„ 

n 

- i ^KejGjF^{co)Ej{co) V G„ /—> 

n V 

with Ej(a>) = ^e‘^x*Rj(~^) + ^~‘^XRj(<^)\ ^^d E*(-a)) = Ej(a)), and have the correlation 


(^i^outj(^)'^i^out,/( ^)) — 2yffiQ ^KgjKglGjGj(^OJ)i^I (m) (^m “1“ 1) “1“ \Xini. ^)l 

1 2-1 2- 
+ 2\'^'^^JXRj(^) ~ 1| (^oj + 1)^;7 + 2^2KejXR,j{-iO) “ 1| flojbjl 

+ e~‘^ ^jKejKelGlGj \2Ke,jXRj{^) “ l] Fm(<^)^K(^)X*Rj((^)(bo,J + 1) 

- e‘^ ^jKejKelGlGj \2KejXR,j{-<^) “ l] Fm(<^)^K^)XRj(-(^)boJ 
+ e‘^ ^jKejKelGlGj ^Ke,lX*Rj{co) - l] Fm(o))Ej(oj)xR,i(a))(no,i + 1) 

- ^[K^iGiGj [2Ke,iXR^i{-io) - 1 ] F„{ix))Ei{io)xR i{-co)no,i 
“1“ 2Ke jKi j\XR,jixY)\ (jloj 1)^2/ 2Kg jKi j\XR,j{, ^)l ^ojbjl 

+ 2Kij^/K;ji^je-‘^GiGj\xRj(a^^^^ + 1) 

- 2Kij ^jKejKejF^GiGj\xR,j(-co)\^F*^(co)E^(co)noj 

+ 2Kij ^jKejKeje‘‘^GiGj\xR,i(co)\^Fm(a>)Ej{a)){no,i + 1 ) 

- 2Kij ^JKeJKeJe~'^GlGj\XR,l(-CO)\^Frn(a))Ej(co)noJ 

■I" 2 ^jKgjKejGlGj\F^(co)\ (cu) ^ ' (yKe,n "I" Ri,n)Gy\XR,niMb)\ {fio,n "I" 1) 

n 

2 "\^Kg jKg iGiGj\Fffi(^oS)\ ,3y'(6U)i3^ (m) ^ “I" ^)l ^o,n • 


( 25 ) 


( 26 ) 
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Merging some terms, we have 



(^) [^U^m(^)l (^m 1) lATmC ^)l 


+ e-‘^ [iKjXRjico) - l] Fl{oj)S;{co)xlj{co){ri,j + 1) 


+ ^jKejKeiGiGj [2/c/^* - l] F^(aj)Ej(aj)xRj((o)(noj + 1) 



(27) 


n 


+ e yjKejKelGiGj \lKjX*R^j{-Oj) - l] F*^((o)Si (oj)XRj(-Oj)noj 


+ e'^ y/KejKeiGiGj [2 kiXr,i(-co) - 1] F,„{i 0 )Sj{(o)x*Rj(-co)noj 


+ 2 ^KejKejGiGj\F^iaj)f-Sj(a))Si(co) ^ KnGl\XR,n{-Oj)\^n„^n 


2 Entanglement 

One of the suffieient eriterion for the bipartite entanglement of any two CV states is Duan’s eriteria [|29ll . Its 
extension for multipartite entanglement of CV states is written directly in terms of the correlation of fluctuation 
of fields, as sums of variances fl^ . Here we derive the variances for any paired CV states to provide an analysis 
for multicolor entanglement of CVs. 

We assume that there are M cavity modes in our optomechanical system involved in operation. Without 
loss of generality, we calculate the sums of variances between two modes aj and a/ by 


=(AVout,;(n>) - AX,„,jiaj))iAX,,,ji-aj) - AX,,,ji-aj)) 

■I" (Ayout.X^) Ayout,/(^))(Ayoutj(~<^) Ayout,/(~^)) • 


(28) 


The sufficient criterion for inseparability between these two modes is Vji < 2. If all bipartition in M CV states 
are entangled, i.e. < 2 for any j I e {1,2, - , M}, genuine multipartite entanglement is then obtained 

mniEii. Note that the way to calculate these variances is not optimal for the criteria for analyzing CV 
entanglement. But the violation of inequality Vji > 2 is sufficient to claim that these two modes are entangled. 
This variance can be smaller if coefficients associated with other modes are properly chosen [ini30ll. 

To provide a fundamental limitation for obtainable entanglement we consider the case all cavity modes have 
identical parameters, i.e. Gy = G, k^j = k^, Kij = ku kj = k and A, = A. We take A ~ 0 and define 6+ = f2,„o ± co 
and Assuming that |d| » ymo, then we have 


2 . Xmi<^) ~--t a.ndxm(-oj) ~ “i;; 
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Using the above approximation we obtain the correlation 




+ ?7- 


(1 - cos20)(2no + 1) 

sin 19 


(Al^outMjAyout,/(-^i^)) “2(2^0 + 1)^;7 + ^-^1 2^ 


+ ?7- 


+ 5;2)(1 - cos 20) + 

O Qm 


(1 + COS 26)(2no + 1) 
sin 26 


(2ro + 1) 


c2/r-2 . c—2\/i . 

— )(1 + cos 20) 


(200 + 1) 


(29a) 


(29b) 


with T] = kJk and fn 


r, and the variance 


y'f^(m) =2(200 + 1) + 4?7M 


20 


(1 + cos 20)(2no + 1) 


477- 


r2/c-2 r-2\/'i o/3\ Sin 20 

— —0 (0_ + 0+ )(1 + COS 20)-— 


(200 + 1 ) . 


(30) 


Obviously, our Eq. ( [29l ) for a single mode, j = I, agrees with the output noise power density in ll2T]| af¬ 
ter replacing k with k/I and setting ho = 0. Entanglement of bipartite CV states is obtained when |0| > 

is minimal, = 2 - 2r] (i + Jp-”"’ ), 

(mT _V MPmeast?/,!/ 


I' vffmeas _|. I I ^ 2 ^ + 1 ) Whllc thc varlancc is minimal, = 2 - In ‘ ^ 

V 2 Q,„ / I sine I ^ o / ' mm 4nm^V 

pwimeas+ j 

at the optimal frequency 0opt ~ (iUEmeas + In the strong coupling regime, Mf^eas ^ , the 

„o,-;o,a^o T/W « 2 - ^ sm^ 0, is indepcndcnt of the optomechanical coupling strength, G. On the 


4 nn,n,, 


minimal variance, V'^ 

conEary, in the weak coupling regime, MT^^as ^ 
proportional to the square of optomechanical coupling strength, G^. 
At 0 = ±7r/4 the variance becomes 


, the minimal variance, V 

■ “ T 


(M) 

min 


rjQmMVrti 


2QinMA}] 


- sin^ 0, is 


vJfV) 


2 + t]M 


rmeas.*^m 


(2no + 1) • 


(31) 


The bipartite entanglement of CV states requires |0| > Mfn,eas/2 -i- 0.^n,„IQ,n. At the optimal frequency |0opt| = 


AfEmeas + 40 


mo^i we obtain the minimal variance, 


'?rn: 


Mfa.ea.+4n.aog 


(200 + 1) > (2 - 
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